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METHOD OF ESTIMATING THE STICK-FIXED LONGITUDINAL 
STABILITY OF WING-FUSELAGE CONFIGURATIONS 
HAVING UNSWEPT OR SWEPT WINGS 
By Milton D. McLaughlin 


SUMMARY 


A method. is given for calculating the stick- fixed longitudinal 
stability of a wing-fuselage configuration at subcritical Mach numbers. 

The method applies to unswept- and swept-wing configurations. No effort 
■ has been made to include the effect of the horizontal tail but this 

effect could "easily be calculated if the downwash values were known. In 
the development of the method experimental results were necessary to 
obtain the normal loading on the part of the fuselage in the vicinity of 
the wing- fuselage . juncture. Experimental curves are presented which are 
used to predict the loading, on this portion of the. fuselage for -other 
configurations. 

The geometrical characteristics and the calculated contributions to 
lift. and pitching moment of the principal parts of the wing-fuselage con- 
figurations for 23 configurations are presented. A comparison of the 
neutral points predicted by the method with the neutral points obtained 
from experiment for the 23 wing-fuselage configurations is also presented. 


INTRODUCTION 


Methods for computing the stick-fixed longitudinal stability of 
airplanes having small sweep angles (below 20°) and high aspect ratios 
have been developed to a point where a good degree of accuracy may be 
obtained with their use. In methods such as that given in reference 1 
the agreement between the calculated and measured neutral points has 
generally been ±1.5 percent of the mean aerodynamic chord or within the 
accuracy with which neutral points have been evaluated from flight data. 
These methods do not apply, however, to configurations having large 
sweep angles. = : . .. 
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The purpose 'of the present paper Is to. derive a method which per- 
mits calculation to within a satisfactory degree of the stick-fixed 
static longitudinal stability of configurations having either unswept 
wings or wings with large sweep angles. This method applies to a wing- 
fuselage configuration. No effort is made to include the effect of the 
horizontal tail hut this effect could easily he calculated if the down- 
wash values were known. Since no generalized charts exist for estimating 
downwash behind swept wings , the downwash values used in the method are 
best obtained from experimental values on a similar wing- fuselage con- 
figuration. Theoretical procedures are used where possible to develop 
this method for calculating the aerodynamic forces and moments; empirical 
methods based on experimental results are used for those parts of the 
configuration for which no theory is available. • 

The stick-fixed neutral points of 23 wing- fuselage configurations 
have been calculated by the present method and a comparison with experi- 
ment is presented. The 23 wing- fuselage configurations together with 
pertinent information- are presented in table I. Results of tests of the 
first 10 configurations provided experimental data which were used to 
obtain the variation of the normal loading on the fuselage in the vicin- 
ity of the wing-fuselage juncture. These experimental curves are used 
to predict the loading for similar configurations. A sample calculation 
of the neutral point of a wing- fuselage configuration is included in an 
appendix . 


SYMBOLS 


b wing span perpendicular to plane of symmetry, feet 

c wing chord parallel to plane of symmetry, feet 

"c mean aerodynamic chord, feet 

c av average wing chord ^S^b), feet 

x longitudinal coordinate measured from quarter chord of mean 

aerodynamic chord, feet 

x D longitudinal distance measured from quarter chord of mean 

aerodynamic chord to duct inlet, feet (positive forward) 


x f longitudinal distance measured from quarter chord of mean 

aerodynamic chord to. a given fuselage station, feet 
(positive forward) 

x f) longitudinal distance measured from fuselage nose to a given 

' fuselage station, feet 
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H 

x i 

y 

Y 0 

y 

Yf 

s 

11 

n 

cp 


X 

A 

M 

V 


q 

a 

e 

P 


longitudinal distance measured from quarter chord of mean 
aerodynamic chord to leading edge of intersection chord, 
feet (positive forward) 

» 

longitudinal distance measured from quarter chord of mean 
aerodynamic chord to center of pressure of carry- through 
section (positive forward) 

lateral coordinate measured from plane of symmetry, feet 

lateral distance to aerodynamic center of external wing 
measured from plane of symmetry, feet 

average fuselage "width at wing root, feet 

finite fuselage width, feet 

area, square feet 

dimensionless lateral coordinate 

spanwise center-of -pres sure location on wing, fraction of 
semispan fefj 


aspect ratio b 


, , . /Tip chord \ 

taper ratio ^ oot choI . a J 

sweep angle, degrees 
Mach number 

velocity, feet per second 

average air velocity at minimum cross section of duct inlet, 
feet per second 

free-stream dynamic pressure, pounds per square foot 
angle of attack, degrees 
downwash angle, degrees 

angle of local air flow relative to x-axis, degrees 
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C Z 

% 

c n 


m 


“a 


C L, 


'a 


'“a 




a„ 




c L - 


a i 


C La T 


factor used in determining dp/ da 

/Local lift\ 


local lift coefficient 


^Local liftj 


local lift coefficient due to additional lift on wing 

^Local normal force j 

/Lift\ 


local normal-force coefficient 


lift coefficient — ) 

\&u ) 

pitching-moment coefficient ( Fltdll °g g oment ) 

variation of local lift coefficient with angle of attack, per 
degree ^dc ^ j da^ 

variation of lift coefficient with angle of attack, per . 
degree ^dCk/da^ 

variation of pitching-moment coefficient with angle of attack, 
per degree |dC m /daj 

variation of center-section (fig, l) lift coefficient with 
angle of attack, per degree (based on wing center- section 
area Sj) 

variation .of external -wing (fig. l) lift coefficient with 

angle of attack, per. degree (based on external -wing area S 0 ) 

variation of external -wing (fig. l) lift coefficient with 
angle of attack, per degree (based on total wing area Stf) 

variation of carry- through- section (fig. l) lift coefficient 
with angle of attack, per degree (based on wing carry- 
through area Ycj_) 

variation of carry- through- section (fig. l) lift coefficient 
with angle of attack, per degree (based on total wing 
area Sw) 


Subscripts: 


o 


external wing 
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W 

WF 

a 

i 

I 

F 

f 

MAC 

D 

EE 

ac 

c/4 

TE 

l 

max 

cp 

av 

exp 


total wing (including center section) 

complete configuration 

center section of total wing (fig. l) 

wing- fuselage intersection 

wing carry- through section (fig. l) 

fuselage (fore and aft) 

finite fuselage section 

mean aerodynamic chord 

intake duct 

leading edge of wing 

aerodynamic center 

quarter chord 

trailing edge of wing 

fuselage length 

maximum 

center of pressure 

average 

experimental 

METHOD OF ANALYSIS 


The stick-fixed neutral point of a wing- fuselage configuration is 
defined as the center -of-gravity location for which the slope of the 
curve of airplane pitching-moment coefficient against lift coefficient 
dC m /dCL is zero. In order to estimate the neutral point of a configura- 
tion it is necessary to determine the additional loading and pitching 
moment of the configuration. In order to determine these quantities 
the configuration is separated into its principal parts and the addi- 
tional loading of these parts is calculated. The resultant forces and 
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pitching moments contributed by each part are presented in coefficient 
form as values of Ci^ and C^, respectively. Prom the summation of 


these coefficients the. distance of the neutral point of the wing- fuselage 
configuration from -the reference axis, expressed as a fraction of the 
mean aerodynamic chord of the wing may be found as the value of 


dC m \ 


where 



( Cna )wr 


In the present method, the configuration is separated into the 
following three parts: the external wing, the fuselage fore and aft 

of the wing-fuselage juncture, and the wing carry- through section. 

(See fig. 1.) Hereinafter, the part of the fuselage fore and aft of 
the wing-fuselage juncture is referred to as the "fuselage." The carry- 
through section is considered rectangular in shape; its length is equal 
to the length of the wing-fuselage intersection chord and its width is 
the average width of the carry- through section. It is important to note, 
however, that whenever the total wing area is considered in the method 
the center-section area is the area indicated by the dotted lines in 
figure 1. The lift and moment coefficients for the principal parts are 
calculated in terms of the total wing area. The moment reference point' 
is taken as the quarter-chord point of the mean aerodynamic chord of the 
total wing. 

The present approach to the analysis of longitudinal stability dif- 
fers somewhat from the approach presented in reference' 1. In refer- 
ence 1 the lift of the wing-fuselage combination is taken as the lift 
of an isolated wing of corresponding dimensions. Fuselage lift, there- 
fore, must be taken as. equal to. the lift on the section of an isolated 
wing which is covered by the fuselage. Fuselage moment is obtained from 
an integration of the moments contributed by sections along the axis of 
the fuselage (equation (3*7), reference 2). The loads on each section- 
are computed on the assumption that the flow is two-dimensional in planes 
perpendicular to the fuselage axis. This assumption is less accurate in 
the neighborhood of the wing and _ doe s not allow calculation of the loads 
on the wing carry- through section. For airplanes with unswept wings, 
however, conditions on the fuselage near the wing do not greatly influ- 
ence the pitching moment of the fuselage about the neutral point of the 
airplane and, as a result, the method is adequate' for such cases. 

For swept-wing configurations the loading on the fuselage and wing 
carry- through section. is much the same as for unswept configurations. 
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The wing carry-through section for swept wings, however, is located 
farther from, the neutral point of the configuration and the longer . 
moment arm results in a large pitching-moment contribution. The lift 
contribution of the wing carry- through section, therefore, must be 
determined as accurately as possible in order to determine satisfacto-. 
rily the pitching-moment contribution. Thus the assumption that the 
fuselage lift can be taken as equal to the lift of the center section 
of an isolated wing no longer applies and the lift and moment contri- 
butions of both the fuselage and the wing carry-through sections must 
be obtained. In general, the method used herein is to estimate the load 
distribution on ‘the fuselage and the load on the external wing by theo- 
retical methods, whereas for the carry- through section, it is necessary 
to develop curves based on available experimental data for predicting 
the loading and the aerodynamic center. The theoretical values which 
were used in calculating the wing stability parameters were obtained 
from reference 3« Reference 3 was used because it presented a readily 
available uniform source of information and showed fairly good agreement 
when checked with experiment. The method to be used in obtaining the 
wing stability parameters, however, is at the discretion of the designer. 

Although the method is derived for incompressible flow, it may be 
applied for any case of sub critical Mach number by use of the generalized 
Prandtl-Glauert law described in reference 4. The effects of Mach num- 
ber on the neutral-point location at subcritical Mach numbers are gen- 
erally small. This method is intended to apply at low lift coefficients 
where the drag effects are least and where the slopes .of the lift curve 
and pitching-moment curve are linear. 


CALCULATION OP STABILITY 


The method of obtaining the values of Ci^ and C^ of the indi- 
vidual parts, as well as of the complete configuration, is now presented. 

,A step-by-step procedure is utilized to facilitate the computations and 
the sample calculations to illustrate the method are given in the 
appendix . 


External Wing 

Calculation of Cj^.- The contribution of the external wing to the 
lift-curve slope Cl^ may be obtained in the following manner: 

(l) Obtain the value of Cl for the total wing from reference 3, 
figure k. 
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(2) Calculate the quantity Cr %. (The value of 3/ K of refer- 

a W 

ence 3 as used in the present paper is considered to be unity.) 

( 3 ) Obtain from reference 3> figure 3^ the spanwise values of 
c 2 a cy / CL c av for the center section of the total wing. Each value of 

c Z a c / A'L c av is equal to the value c l a C / < ''i J cfy Cav ’‘ 


(4) Multiply each quantity c-j c/C^ c av by the product of the 

mean geometric chord c av and the liftr curve slope for the total wing 
C T to obtain, the quantity c 7 c. 

Lq W a 

(5) Plot the values of c^ c against their spanwise location for, 
the center section of the wing. 


(6) Integrate the value C7 c spaqwise over the wing center sec- 

CL 


tion. The integral is equal to the quantity Ct S ff . 


(7) Obtain the quantity Cx,t S D for the external wing by the 

“o 


formula 


°i\ 3 ° = ■ Cl % s ° 


(8) Obtain the lift- curve slope of the external wing based on 
total wing area by the relation 


Cl 


'cl 


Cl 'o v S ° 

% 


Calculation of Cm^.- The contribution of the external wing to 
pitching-moment-curve slope Cm^ is calculated as follows: 


(l) Obtain from reference 3* figure 5, the spanwise center of 

pressure of the lift of the total wing t) as a fraction of the 

cp 

semispan. - * — ‘ 


MCA EM L5U23 


9 




Y/2 


( 2 ) Perform the spanwise integration indicated in the function 


c Z a cy dy 




Y 72 " 


over the semispan of the wing center section to obtain 


c 7 c dy 
‘a 


the spanwise center of pressure of the wing center section. In most 
cases the values of c^c are nearly constant over the wing center sec- 
tion and the center of pressure is located at one-half the center- 

1 Y 

section semispan, that is y a = o o‘ 

cp t - 


(3) Obtain the spanwise center of pressure of the external wing 
from the relation 


cp 


cp 


Cr, SuYu - Ct S_y_. 
ij cfy w w cp "“otjj a °cp 

Cr Su - Cl S/t 

. -aw % 


This relation is the spanwise center of pressure of the total wing 
corrected for the removal of the lift on the center section. 

(4) Obtain the contribution of the external wing to the pitching- 
moment slope. Cm by the relation 




te* A c/4 C L 

' °o 


If the aerodynamic center of the external wing does not lie on the 
quarter chord (as might be the case if other than the lifting- line 
theory of reference 3 were used), the distance between the quarter- 
chord line and the aerodynamic center must be considered in the preceding 
relation. The pitching-moment contribution of the external wing is then 
obtained by the relation 



/ y MA0 - y °cp 



c 


tan A c / k + 
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Fuselage (Fore and Aft of the Intersection Chords) 
Calculation of Cl^.- The contribution of the fuselage to lift- 


curve slope C; 




is determined as follows: 


(l) Obtain the contribution to lift- curve slope of the fore and 
aft fuselage sections by a graphical integration of the following 
formula which was derived from formula (3*3) ^ reference 2: 




L 


LB d 





where Yf is the width of fuselage at any station. The fuselage plan 
form is divided into a finite number of sections as shown by figure 1. 

For each of these sections the parameter Yf^ ~ is calculated and 

plotted against its fuselage longitudinal station x.p t . The measured 

\2 Y ? ai) 


slope 


dxj 


at each of these stations is proportional to the 


slope of the sectional lift coefficient cj (per radian) of that 

Of 

station. 

For the fuselage stations ahead of the wing values of dp/da are 
obtained from figure 2 and from the following formula: 


= 1 + k 57.3 

da A 


which is the variation in local air flow with angle of attack for the 
fuselage alone plus the additional variation caused by the presence of 
the wing. Figure 2 was derived from theoretical data obtained from 
reference 5« The parameter k is noted to' be expressed as a function 
of the distance forward of the intersection quarter chord c^y^ which 

should not be confused with the quarter chord of the mean aerodynamic 
chord "c/b. For the fuselage sections behind.' the wing, values of dp/da 
are obtained from the variation of downwash with angle of attack de/dtt 
by the formula 


dg _ 1 _ d£ 
da da 
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As no generalized charts exist for estimating downwash ■‘behind swept 
wings, experimental values of de/da at the location of the tail 
should he used when possible. The values of de/da are assumed to 
vary linearly from, zero at the wing trailing edge to the determined 
value at the aft end of the fuselage. The exact value of downwash at 
the aft end of the fuselage is of little importance in determining 
fuselage lift or moments inasmuch as the lift and moment contributions 
of the rear section of the fuselage are very small. 


Calculation of Cj^.- The contribution of the fuselage to > 

pitching-moment-curve slope Cm^ includes the following steps: 

(l) Obtain the contribution of the fore and aft fuselage sections 
to pitching-moment- curve slope by a graphical integration of the fol- 
lowing formula which was derived from formula (3*7 ) > reference 2: 


C 


- 1 

% " 57- 3Sjj5; 



Xf 


dxf i 



where Xf is the longitudinal distance measured from the quarter 
chord of the mean aerodynamic chord to a finite fuselage station and 
x~, is the longitudinal distance measured from the fuselage nose to a 
finite fuselage section. 

(2) Obtain, when appropriate, the additional contribution to 
pitching-moment-curve slope of a fuselage air-intake duct from the 
following formula: 


1 V D 2SdX D 

“op ~ 57-3 V s 

The effect on pitching moment of the change in momentum of the air 
entering the Id take duct obtained from this formula is illustrated in 
figure 3* The formula was derived from a similar unpublished formula . 


Wing' Carry-Through Section 

Calculation of Cp^.- The contribution to lift- curve slope of the 
wing carry- thro ugh section Cp^ is calculated as follows: 
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(1) Obtain the contribution to lift- curve slope of the wing carry - 
through section from figure 4 ( 3 ) by using the effective aspect ratio of 
the carry- through section. The lift-curve slope of the carry- through 
section is based on that area (intersection chord Cj_ times average 
fuselage width at wing root Y) and is presented in terms of the 
section-lift-curve slope of the wing at the intersection chord 

as Cl ./c*^. The value c 2 ^ may be obtained from step ( 5 ) of the 

outline of the contribution of the external wing to lift-curve 
slope C]\ » 

no 

(2) Calculate the contribution to lift-curve slope of the wing 
carry -through section by the relation 



C L « 


<*1 


ta i 


c z 


°i 


Calculation of Cm^.- The contribution to pitching-moment- curve 
slope of the wing carry-through section Cm^ is obtained as follows: 


(l) Obtain the aerodynamic center of the carry- through section 
from figure 4(b) , which is a plot of the aerodynamic center of the 
carry- through section against sweep angle of the quarter-chord line of 
the wing. The aerodynamic- center location is expressed in intersection 
chords as the longitudinal distance from the leading edge of the inter- 


section chord 



(2) Calculate the moment arm of the aerodynamic center of the 
carry- through section about the reference point by the relation 



( 3 ) Calculate the contribution to pitching moment of the carry' 
through section by the relation 
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Complete Configuration 

Calculation of Cr Obtain C r as a sum of the lift-curve 

_^F 

slopes of the principal parts 

Cl owf ” ' Cl cc 0 + Cl c^. + CL aj 

Calculation of Cm Determine C™ by summing up the 

5?F ^WF 

pitching-moment- curve slopes of the principal parts 

Cm a WF = C % 0 + C “c^ + Cm aj 

Calculation of neutral -point location .- Obtain the neutral-point 
location with respect to the reference axis, expressed as a fraction of 
the mean aerodynamic chord of the wing by the relation 

Cmq WF _ fWF 
C L-* c 

DISCUSSION OF METHOD 


The values of wing lift-curve slope for the total wing were 
obtained from reference 3 and are based on the lifting- line method of 
calculation. Values of wing lift-curve slope were also calculated by 
the method of reference 6 which is essentially a correction of unswept- 
wing data for sweep effects. This method has the advantage of allowing 
the use of values of lift-curve slope for unswept wings based on 
lifting- surface theory. Of the two methods the lift-curve-slope values 
from reference 3 showed the best over-all agreement with the experimental 
values for the particular configurations checked in this analysis. 

In order to determine the lift-curve slope of the external wing 
the sectional lift is assumed to be the same as the sectional lift for 
a corresponding section on an isolated wing. Upflow around the body 
might be expected to increase the lift on the sections near the wing 
root. Experimental investigations of this effect are rather scarce, 
however, the data of reference 7 indicate that, if this expected 
increase does occur, it is very small. Also, flight measurements of 
the spanwise center of load of a wing in combination with a body have 
shown good agreement with the values predicted on the assumption that 
the sectional loading external to the body is the same as that on an 
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isolated wing. The lift-curve slope of the external wing therefore, is 
assumed to be equal to the integration of the spanwise loading on the 

Sq 

external wing in terms of the total wing area C^r — . The 

a o % 

aerodynamic -center location for the external wing is then the 
aerodynamic -center location of the total wing corrected for the removal 
of the center section. 

The present method for determining the values of and of 

the fuselage is based on Multhopp's method in reference 2. Good agree- 
ment is found to exist between Multhopp's method and experiment on the 
fore and aft sections of the fuselage as is shown in figure 5- The 
experimental fuselage sectional loading (fig. 5) was obtained from 
unpublished pressure-distribution measurements on a wing- fuselage con- 
figuration (configuration 10 in table I). For the fuselage section 
adjacent to the wing large disagreement can be seen between theory and 
experiment. The large difference in loading on the sections immediate 
to the wing was caused by the wing and in the present method this dif- 
ference in loading is considered to be part of the wing carry-through 
loading. 


Multhopp's method does not apply for the wing carry- through section. 
An effort was made, therefore, to obtain from experiment the loading of 
the carry- through section and its possible variations with the wings of 
various sweep, angles and aspect ratios. These experimental curves 
could then be used to predict the loading on similar carry-through sec- 
tions for other configurations. A limited amount of experimental 
pressure-distribution data was available on loadings of the carry- 
through section; however, a larger amount of data was available on the 
over-all characteristics of wing- fuselage configurations. From the 
wing-fuselage data and the isolated-wing data the experimental values 
for the wing-fuselage and total wing were obtained. By using the pres- 
ent method, the contribution to stability of the fore and aft fuselage 
sections was calculated and the effect of the wing center section was 
calculated and subtracted from the experimentally measured contribution 
of the total wing. The contributions to stability of the wing carry- 
through section were then found as 



K>Lp + M theory 





theory 
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Experimental values were used in calculating the contributions to 
lift and pitching moment of the external wing for the following reasons: 
The contribution to lift of the external wing is large when compared 
with the contribution of the wing carry- through section; consequently, 
a small error introduced in the wing contribution to lift results in a 
much larger error when expressed in terms of the contribution of the 
carry- through section as can be seen from the formula. It was therefore 
believed that, by use of experimental values of wing lift and 
aerodynamic-center location in the preceding formulas, a more accurate 
contribution of the carry through- se'ction lift and, hence, to moment 
could be obtained. Table II presents the experimental value of 


and 

and 


CL. 

Cm aj* 


which were used in calculating the values of 



A plot of values of 



(based on carry- through area) with 


sweep angle of the quarter- chord line is shown in figure 6. This 

lift-curve slope of the carry-through section appears to be fairly con- 
sistent about a value of 0.055 per degree for the configurations checked. 
A somewhat better correlation with experiment was obtained, however, by 
another method of analysis. A plot of Cti fci (the lift- curve 

®c/ °i 


slope of the carry- through section in terms of the intersection-chord 
lift-curve slope) against l/c± (fig. 4(a)) indicates that the loading 
of the carry- through section varies with the loading of the adjacent 
wing and also varies inversely with effective aspect ratio of the carry- 
through section. A check point obtained. from actual integration of the 
experimental loading of the carry- through section presented in figure 5 
iB shown in figures 4 and 6. 


An examination of figure 4(a) reveals that, as the effective 
aspect ratio of the carry- through section approaches zero, the value of 
the ratio Cr > fci exceeds the value of unity. This fact may be 
a I/ °i . 

explained by an examination of figure 5» The loading of the carry- 
through section is found to include part of the adjacent-fuselage 
loading fore and aft of the carry- through section and, since this 
loading is based on the carry- through area, the wing- induced fuselage 
loading may possibly exceed the adjacent intersection-chord loading. 

When the effective aspect ratio of the carry-through section approaches 
zero, however, the value of the ratio Ct t /ci must approach unity. 

a I/ «i y 

In the absence of experimental data in the region below — = 0.325 

c i 

(fig. 4(a)) no attempt has been made to fair the curve to zero effective 
aspect ratio. 
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The chordwise aerodynamic-center location of the carry- through 
section is shown in figure 4(h) as a plot of 


/* i\ _ fa) 

l c i/LE \ c iy 

from the leading edge of the intersection chord) against 


(measured 

(&c 

Shown 

also in figure 4(b) is the chordwise variation with sweep of the 
aerodynamic -center location for the center chord as given by the flat=“ 
plate first-approximation theory (reference 8). The experimental values 
show approximately the same variation with sweep as the flat-plate theory 
except at high- sweep angles. The experimental values also tended to 
show some effects of the .Reynolds number. For a constant sweep angle 
a higher Reynolds number resulted in a more rearward location of the 
computed aerodynamic center. Because of the lack of additional experi- 
mental points, only one curve has been faired in figure 4(b) . The 
check point obtained from actual .integration of the loading of the 
carry- through section (fig. 5) is also shown in figure 4(b) and shows - 
good agreement. 


The configurations used in obtaining the loading and aerodynamic 
center of the carry- through section shown in figures 4 and 6 and the 
check configuration are listed as the first 10 configurations in tables I 
and II. For these 10 configurations sufficient experimental data were 
presented to obtain the loading on the carry- through section. The 
static longitudinal stability of the 23 wing- fuselage configurations 
was calculated by using the curveB faired through the experimental data 
in figure 4 and by the present method. 


CORRELATION OF RESULTS 


The method has been used to determine the static longitudinal 
stability of the 23 wing-fuselage configurations listed In table I and 
the results are presented and compared with the experimental results in 
table II and figures 7 to 9* . A comparison between the experimental and 
calculated values of lift-curve slopes is presented in figure 7* The 
calculated values show some disagreement with experiment, however,, the 
disagreement was found to result mostly from low values of CLq^ 

(theoretical) for wings with approximately 45° sweep angle. 

A comparison between the experimental and the calculated values of 
pitching-moment-curve slopes for the configurations is presented in 
figure 8. The results of the calculations show reasonable agreement 
with the experimental values. The disagreement which exists, however, 
can be traced primarily to the disagreement in lift-curve slopes and 
to the limitations in estimating the aerodynamic- center location for 
some of the lift contributing parts. Calculation of the aerodynamic- 
center location of the wing carry-through section shows that there are 



MCA EM L5U23 


17 


probable Reynolds number effects which have not been considered in the 
method. 

A comparison between the experimental and calculated neutral points 
for the various configurations is presented in figure 9* The agreement, 
which is generally better than ±0.04, is considered to be good, espe-' 
cially for such a large variety of wing plan forms. The accuracy with 
which the aerodynamic center and the lift- curve slope of the wing can 
be determined to a large extent determines the accuracy of the neutral- 
point location for a wing-fuselage configuration. 


CONCLUDING REMARKS 


A method is given for calculating the stick-fixed longitudinal 
stability of a wing-fuselage configuration with unswept or swept wings 
at subcritical Mach numbers. The stability parameters estimated by the 
method show reasonable agreement with the experimental values for the 
23 configurations used in the comparison. Calculated values of lift- 
curve slope show some disagreement with experiment at sweep angles close 
to 45°. This disagreement is found to result mostly from low theoreti- 
cal values of wing lift-curve slope. The effect of the low values of 
wing lift-curve slope is also apparent in the calculated values of the 
slope of the pitching-moment curve. There are noticeable Reynolds num- 
ber effects on the aerodynamic- center location for the wing and wing 
carry- through section. As no definite trend of these effects was 
established, the effects are not considered in the method. To a large 
extent the accuracy with which the wing stability parameters can be 
obtained determines the accuracy to which the configuration stability 
parameters can be calculated. 

For the wing carry-through section no theory was available for 
predicting the loading; therefore, the loading was obtained from curves 
based on experiment. The experimental lift-curve slope of the wing 
carry-through section calculated in terms of the carry-through area 
appears to.be consistent about a value of 0.055 per degree. A better 
correlation with experiment, however, indicates that the carry-through- 
section loading varies with the loading on the adjacent wing section 
and also varies inversely with the effective aspect ratio ' of the carry- 
through section. 

The method has been used to calculate the neutral points of 23 con- 
figurations. The agreement between experimental and calculated values 
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has generally been better than ±0.04 and is considered to be good 
especially for such a large variety of-ving plan forms. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPENDIX 

SAMPLE CALCULATIONS 




The contributions to lift-curve and pitching-moment-curve slopes 
and the aerodynamic-center location were calculated for wing-fuselage 
configuration h (see tables I and H) by the method presented in this 
paper. These calculations are presented as an aid to further under- 
standing the method. The steps are numbered to correspond with the 
steps presented in the section "Calculation of Stability." Two views 
of the model are shown in figure 10 which was obtained from reference 9 • 
Some total wing parameters used in the example are as follows: 


= 3-5 

Y = 3-6 feet 

= 27 feet 

y MAC = 5 ^ feet 

- 8 .6k feet 

\ = 0.25 

=7.72 feet 

A /, = 60 . 8 ° 

c/h 

= 208.3 square feet 


External Wing 

Calculation of C^ .- Steps used in calculating are: 

(1) From reference 3> figure k, Cl^ = 0.0^12. 

t 

( 2 ) Cr _.S w = (0.01H2) (208.3) = 8.60 square feet 

(3) From reference 3> figure 3: 

For T) = 0 


C 7 C 

4 a 

c L c av 


= 1.13 = 


c Zq c 

'Lo^cav 
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For T) = 0.13 


C 7 C 


c Z r C 


a =1.15= ' ta 


c L c av 


Ct ’ C av . 


C Z, C 


(4) Ct C = -rq — Ct C av 

^ C Lo w c av Lq W av 

For tj = 0 


Ct C = (1.13) (0.0412) (7.72) = 0.359 feet 
4 a 


For r] = 0.13 

c] c = (1.15) (0.0412) (7.72) = 0.365 feet 


a 


(5) See figure 11 for plot of local wing loading c 2 a c 
wing spanwise station. 


(6) 0^ = (c Ia =) av Y 


= °‘3?9 * . °:3 . S 3.6 = 1.302 square feet 


(7) C LV S o = C L nri s W - . c L a S (T 


(8 ) Cl 


= 8.60 - 1.302 = 7-298 square feet 
Ct i S_ 


Oo % 

- 2^3 = O' 0349 

Calculation of C™ .- Steps for calculating C™ are: 



(l) From reference 3> figure 5(t>) > ti = 0.449. 


y w = | “ 0.449(13.5) = 6.06 feet 

cp = ^ 


against 
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(2) * I I 


cp 


= | ( 1 . 8 ) 


= 0.9 foot 


(3) y 0 = 
°cp 


~ 8-y-. 

\ w w cp L a g q a cp 
CL^Sw - Ci^Scr 


l= (8.6) (6.06) - (1.302) (0.9) 
8.6 -1.302 


=6.99 feet 


W C 


y MAC - y 0 


C £ 


tan A, 


/^ C L r 


c “ / “Oq 

= (1.789) (0.0349) 


= -0.0115 


Fuselage 


Calculation of C 


^05- 


and C 




Calculation of the fuselage 


sectional lift is shown in table III and plotted in figure 12. The 
fuselage contribution to lift and pitching moment obtained from * an 
integration of the fuselage sectional lift in figure 12 is found to be 


C^ = 0.00366 


Wing Carry-Through Section 

Some wing carry- through- section parameters are as follows: 

Y 

~ = 0.32 

c i 

Sj = Yc^ = (3.6) (11.25) = 40.5 square feet 

x i = 9.33 
c I ai = °*°325 . 



22 


NACA EM L5U23 


Calculation of 


C t.^ . - Steps for obtaining Cy 


follow: 


( l) From figure 4( a ) , 


C L 


c l 


°i 


1.42 


(2) C L 

a I 




C Z 

«i 


fl 

Sw 


lin 5 

= 1.42(0.0325) 20573 = 0.00896 

Calculation of C^ Steps for calculating C,^ are: 


(l) From figure 4(b) , =0.68 

V Ci 'tE V Ci 'ac 



= (0.15) (O.OO896) = 0.00175 


Complete Configuration 

Calculation of . - Values of lift-curve slopes for the 

complete configuration are: 

Estimated value by present method 
CLc£ WF * CLa o + CLa F + CLa I 


= 0.0349 + 0.00177 + O.OO896 = 0.0456 
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Experimental value (table H) 

- °- oh i 

Calculation of The pitching-moment- curve elopes for the 


complete configuration are: 

Estimated value by present method 

Cm = Cm + C m + Cm 
^WF “Oq “Op “Oj 

= -0.0115 + 0.00366 + 0.00175 = -0.00609 

Experimental value (table II) 

c m C r wr " 

Calculation of neutral -point location .- The neutral-point location 
is calculated to be: 

Estimated value by present method 

•Stf *HP 


" La WF 


-0.00609 n 1 -ii. 

oToSsS * “°- 13 ^ 


Experimental value (table II) 


Umq WF _ X WF 
Cl o wf c 


-0.06l 

0.047 


= -0.13 
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TABU I - 


DIMENSIONS AND INFORMATION OK CONFIGURATIONS 


Configuration 

Configu- 

ration 

number 

Symbol 

NACA 

airfoil 

A c/4 

(deg) 

\ 

A 

b 

(ft) 

(«> 

1 

(ft) 

(«) 

Sw 

(sq ft) 

Reynolds 

number 

Reference 


F 

1 

O 

65-130 

0 

0.50 

6 

7 

1.202 

1.14 

1,42 

8.125 

1.94 X 10 6 

10 

A 


2 

□ 

64-210 

45 

.383 

5.1 

12.46 

2.6 

1.4 

3.28 

30.35 

6.0 

11 

i 

A 


3 

O 

65-110 

40 

.42 

2.5 

5-83 

2 . *7 

1 . 1*7 

2.91 

13.64 

2.87 

12 



4 

A 

a &lA006 

60.8 

*25 

3-5 

27 

8.639 

3.6 

11.25 

208.3 

8.00 

9 



5 

V 

64 1 -102 

50 

.625 

2.88 

9.41 

3.33 

1.4 

3-8 

30.75 

6.80 

13 

A 


6 

> 

0012 

45 

1 

2.63 

3.04 

1.18 

0.5 

1.18 

3-50 

l.4o 

14 

A 


7 

< 

64^-112 

40 

.625 

4.01 

11.35 

2.89 

1.4 

3 . 3 * 

32.24 

3.04 

15 

A 


8 

17 

6^-112 

40 

.625 

4.01 

11.35 

2.89 

1.4 

3 . 3 * 

32.24 

6.84 

15 and 16 

A 


9 


64 1 -112 

4o 

.625 

*.01 

11,35 

2.89 

1.4 

3 . 3 * 

32.2* 

8.09 

15 and 17 

■ 


10 

E\ 

a 65A006 

45 

.6 

4 

2 

0.51 

0.267 

0.59 

1.00 

1.73 

Unpublished 

X 

Y 

11 

A 

65-110 

0 

.445 

6.8 

7.68 

1.181 

1.14 

1.42 

8.67 

1.27 

10 

X 

h 

12 

D 

65-110 

15 

.46i 

6.4 

7.29 
i 

1.20 

1,10 

1 

1.43 

8.40 

1.292 

i 

10 


a Airfoil par&llel to plane of syimetry (airfoils not- footnoted perpendicular to 0. 25-chord line) . 
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TABU I - Concluded 


DIMENSIONS ARB UffOEHATIOT ON CCNFIGORAIIOifS - Concluded 



“Airfoil perpendicular to 0- 50-chord line, ( airfoils ncrt-footnoted perpendicular to 0 -25-chord line). 


c Root chord; tip chord. 
^Republic airfoil section. 























































































Configu- 

ration 

number 

I 

Symbol 1 

1 

O 

a 

□ 

3 

O 

, 4 

A 

5 

V 

6 

> 

T 

< 

0 

7 

9 


10 

bw 

n 

4 

12 

V 

13 


14 

k- 

30 

a m 

16 

G 

17 

a 

10 

D 

19' 

a 

BO 

Q 

21 

c 

22 

0 

23 

D 


0.074 o.sBa 
.o6e -.0733 

.o 46 .007 

‘ ,p 4 a -.13 
.046 -. 00 G 5 

.046 .054 

.060 .000 

.060 -.003 

.661 .0143 


0.0597 
.0501 
.0330 
- 03 *9 
.0366 

■ 

. 01*66 

.0466 

. 01*66 

.0448 

.0646 

.0616 

.0529 

.0392 

.0503 

.0560 

.053 

.0461 

.o6e0 

.0524 

.0353 

.0062 

.038 


0 

-.0094 

-.oo 4 o 

-.0105 

-,oo 4 i 
-.0005 
-.0046 ' 
-.0046 
-.0046 

-.0057 

0 

-.0037 

-.0068 

-.0000 

0 

-.0007 

-.0020 

-.0046 

-.0033 

-.0056 

-.0068 

-.0046 

.0005 


Calculated 

Calculated 

Cl c* 

Cbq f 

0.0096 

0.0132 

.0022 

.0050 

-0034 

.0045 

.OOlB 

.0037 

,0021 

.0036 

,0022 

' .00£j6 

-0022 

.0043 

~ .0022 

.0043 

' .0022 

' .0043 

.0027 

.0051 

.009 

.0126 

‘ .0078 

.0140 

.0067 

.0140 

.0060 

.0119 

.0038 

.0037 

.0034 

.0046 

.0033 

.0054 

.0030 

.0066 

.0030 

.0048 

.0025 

.0038 

.0022 

.0027 

.0021 

.0018 

.0029 

.oo 4 e 



Esri^i m vDVfi 
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It V 2 


Fuselage 

x f 

Y f 

2 T f 

*f - ^Ci/U 

station 

(ft) 

(ft) 

(eq. ft) 





Iy 2dP 

K 

M 

2 x f da 

da 

(sq ft) 


Fore fuselage 




Aft fuselage 


10 . 

-3.47 

3.50 

19.300 

11' 

-5-47 

■3.34 

17.550 

12 

-7.47 

3.14 

15.500 

13 

-9-47 

2.90 

13.210 

14 * 

H 

i — 1 

i 

2.60 

10.600 

15 

-13-47 

2.24 

7.900 

16 

-15. 47' 

1.82 

5.210 

17 

-17.47 

1.28 

2.580 

IS 

-19.47 

.60 

.566 

19 

-20.47 

0 

0 



25.53 

0 

0 

19.15 

0.710 

0.020 

-4- 

H 

O 

H 

0 

0 

24.53 

.60 

.566 

18.15 

.672 

.021 

1.014 

.575 

.71 

22.53 

1.28 

2.580 

16.15 

•598 ■ 

.022 

1.015 

2.620 

1.26 

20. ; 53 

1.82 

5.210 

14.15 

.524 

.028 

1.019 

5.310 

1.38 

18:53 

2.24 

7.900 

12.15 

' .450 

.033 

1.022 

8.080 

• 1.39 

16.53 

2.60 

10.600 

10.15’ 

.376 

.042 

1.028 

10.900 

1.39 

14.53 

2.90 

13.210 

8.15 

.302 

.055 

1.037 

13.700 

1.39 

12.53 

3-14 

15.500 

6.15 

.228 

.077 

1.052 

16.320. 

.1.25 

10.53 

3.34 

17.550 

4.15 

.189 

.093 

1.063 

18.680 

1.35' 

8.53 

3.50 

19.300 

2.15 

.080 

i 

.245 

1.165 

22.500 

2.55' 
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Figure 1.- Wing-fuselage configuration illustrating separation of the 
configuration into the principal parts.. Some general notations are 
included. ' 
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Distance from. centroid of finite fuselage section to guar ter chord of 

xf - (xci) e /i. 

intersection chord expressed in wing spans, ZLZ 

b 


Figure 2.- Variation of values of n with distance forward from the 
quarter-chord point of the wing-fuselage intersection chord. 

(Reference 5.) Values of dp/da may he obtained for wings of various 

hr 

r i ■ it «.n i t _ -t i ” f CfT 


aspect ratios and liTt-curve slopes ty — \= 1 + k — — 


(57-3). 


V 



Moment reference 
point 




Figure 3 *- Diagram shoving velocity vectors for calculation of addi- 
tional pitching-moment contributions contributed by the air mass on 
the inlet duct. (Unpublished data.) 
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Wing carry- through, section 
-curre elope ®* r ( ^ e S ree 



Check: point 


Sweep angle of quarter chord, A^, deg 


Figure 6.- Wing carry- through, section lift-curve slopes ("based on wing 
carry-through section area Sj) for 10 wing-fuselage configurations 
as obtained by the present method. Symbols refer to configurations 
of tables I and II . Check point refers to actual integration of 
carry- thro ugh loading. 
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Experimental values of lift- curve slope, Cl^, per degree 


Figure 7*- A comparison "between experimental and theoretical values 
of ^ or ^3 wing-fuselage configurations. Symbols refer to 

configurations identified in tables I and II . 




Calculated values of pitcMng- moment- curve slope 


« 


0 


t 



Figure 8.- A comparison between experimental and theoretical values 
of for 23 wing-fuselage configurations. Symbols refer to 

configurations of tables X and II. ^ 
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Figure 11.- Spanwise load distribution of wing semispan for 

configuration 4. 


NACA EM L5U23 



NACA EM L51J23 






Figure 12.- Calculated fuselage lift distribution. 
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